FERMILAB-PUB-05-498-A 



Non-Gaussianity from Broken Symmetries 
Edward W. Kollfl 

Particle Astrophysics Center, Fermi National Accelerator Laboratory, Batavia, Illinois 60510-0500, USA 
and Department of Astronomy and Astrophysics, Enrico Fermi Institute, 
University of Chicago, Chicago, Illinois 60637-1433, USA 



in 
o 
o 

(N 

> 
O 



Antonio Riottcy 

CERN Theory Division CII-1211, Geneva 23, Switzerland 

Alberto Vallinottcjl 

Physics Department, The University of Chicago, Chicago, Illinois 60637-1433, USA 
and Particle Astrophysics Center, Fermi National Accelerator Laboratory, Batavia, Illinois 60510-0500, USA 

(Dated: February 5, 2008) 

Recently we studied inflation models in which the inflaton potential is characterized by an under- 
lying approximate global symmetry. In the first work we pointed out that in such a model curvature 
perturbations are generated after the end of the slow-roll phase of inflation. In this work we develop 
further the observational implications of the model and compute the degree of non-Gaussianity pre- 
dicted in the scenario. We find that the corresponding nonlinearity parameter, /jvl , can be as large 



as 10^ 
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I. INTRODUCTION 

Inflation has become the dominant paradigm to un- 
derstand the initial conditions for the density perturba- 
tions in the early Universe, which are the seeds for the 
Large-Scale Structure (LSS) and for Cosmic Microwave 
Background (CMB) temperature anisotropics 0. In the 
inflationary picture, primordial density and gravity-wave 
fluctuations are created from quantum fluctuations and 
"redshifted" out of the horizon during an early period of 
superluminal expansion of the Universe. 

Despite the simplicity of the inflationary concept, the 
mechanism by which cosmological curvature (adiabatic) 
perturbations are generated is not yet established. In the 
standard slow-roll inflationary scenario associated with 
a single inflaton field, density perturbations are due to 
quantum fluctuations of the inflaton itself. In the cur- 
vaton mechanism 0, the final curvature perturbation C 
is produced from an initial isocurvature mode associated 
with quantum fluctuations of a light scalar (other than 
the inflaton), the curvaton, whose energy density is neg- 
ligible during inflation and which decays much after the 
end of inflation. Recently, we have proposed an alterna- 
tive scenario for the generation of the comoving curvature 
perturbation 

The idea is that if the inflaton sector is character- 
ized by an approximate global symmetry, the associated 
pseudo-Nambu-Goldstone boson(s) may be quantum me- 
chanically excited during inflation. This results in differ- 



* Electronic address: |rocky@fnal.gov|| 
tElectronic address: antonio.riotto@pd.inf riTit| 
■f Electronic address: avallinol 



ent initial conditions in separate horizon volumes for the 
(pre)heating stage after the end of inflation. These differ- 
ent initial conditions will give rise to fluctuations in the 
comoving number densities of the light relativistic states 
produced during the decay process and, ultimately, to 
CMB anisotropics. Other mechanisms for the generation 
of cosmological perturbations have been proposed as well, 
such as the inhomogeneous reheating scenario Q , ghost- 
inflation , and the D-celeration scenario |3| , to mention 
a few. 

Given all these theoretical models, it is very important 
to find ways to discriminate among them based on obser- 
vational evidence. A precise measurement of the spectral 
index of comoving curvature perturbations will pro- 
vide a powerful constraint to slow-roll inflation models 
and to the standard scenario for the generation of cosmo- 
logical perturbations, which predicts — 1| significantly 
below unity. However, alternative mechanisms generally 
also predict a value of riij very close to unity. Thus, even 
a precise measurement of the spectral index will not al- 
low us to discriminate efficiently among them. Also, the 
lack of gravity- wave signals in CMB anisotropics will not 
provide any information about the perturbation genera- 
tion mechanism, since alternative mechanisms predict an 
amplitude of gravity waves far too small to be detectable 
by future experiments aimed at observing the i3-mode of 
the CMB polarization. 

There is, however, a third observable which will prove 
fundamental in providing information about the mech- 
anism chosen by Nature to produce the structures we 
see today. It is the deviation from pure Gaussian statis- 
tics, i.e., the presence of higher-order connected corre- 
lation functions of CMB anisotropics. Since for every 
scenario there exists a well defined prediction for the 
strength of non-Gaussianity and its dependence on the 
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parameters, testing the non-Gaussian level of primordial 
fluctuations is one of the most powerful probes of infla- 
tion 0, and is crucial to discriminate among different, 
but otherwise observationally indistinguishable, mecha- 
nisms. For instance, the single-field slow-roll inflation 
model produces negligible non-Gaussianity, and its dom- 
inant contribution arises from the evolution of the ubiqui- 
tous second-order perturbations after inflation, which is 
potentially detectable with future observations of CMB 
temperature and polarization anisotropies. This effect 
must exist regardless of the specific inflationary model 
considered, setting a lower bound of the non-Gaussianity 
of cosmological perturbations. If evidence for this ubiq- 
uitous non-Gaussianity is not found experimentally, the 
present understanding of the evolution of cosmological 
perturbations will be challenged at a very deep level. 

Motivated by the extreme relevance of pursuing non- 
Gaussianity in the CMB anisotropies, in this paper we 
calculate the level of non-Gaussianity at large scales in 
the CMB anisotropies produced in the class of models 
where the comoving curvature perturbations are gen- 
erated by broken symmetries in the inflaton sector 
Since we will flnd that in this class of models the non- 
Gaussianity is substantial, as a first approximation it 
will be possible to neglect the non-linearity generated 
by gravitational dynamics in the post-inflationary phase 

The paper is organized as follows. In Sec.^]we present 
an overview of the general physical mechanism leading 
to curvature perturbations production during multi-field 
inflaton decay. Sec. IIIII then focuses on the determina- 
tion of the comoving curvature perturbation produced 
when the inflaton field is characterized by a broken U(l) 
symmetry and when the decay process is the so-called 
instant preheating. In Sec. IIVI we carry out the calcu- 
lation of the non-linearity parameter /jvl for the model 
considered. Sec. then generalizes such results to de- 
termine a general formula for the evaluation of the level 
of non-Gaussianity produced whenever the potential is 
characterized by a broken symmetry, and summarizes the 
conclusions. 



II. MODEL OVERVIEW 

As a first step, let's analyze the consequences of as- 
suming an inflationary potential characterized by a global 
broken symmetry for the production of curvature pertur- 
bations during the decay process of the inflaton field. 

In what follows, we consider a multi-field inflationary 
scenario where the scale factor a{t) and the background 
inflaton flelds (pi {t) evolve according to the system of cou- 
pled differential equations 



dV 

3H(j)i + —— ^0, i = l,...,n, (la) 

0(Pi 

8. r^i+^J, (lb) 



where 0= (0i, ^2, f/)™)- 

As in the standard single-field scenario, the multi-field 
scenario is characterized by an initial slow-roll phase dur- 
ing which quantum fluctuations of all the light degrees 
of freedom are excited, stretched to macroscopic scales, 
and frozen in amplitude once their wavelength exceeds 
the Hubble radius. It is important to note that in the 
single-field scenario only fluctuations parallel to the di- 
rection of motion in field space could be generated. As 
was shown in Refs. 0, 0| , these fluctuations would later 
correspond to adiabatic perturbations. In the multi-field 
case, on the other hand, fluctuations both parallel and 
orthogonal to the direction of motion in field space can 
be excited, the latter corresponding to isocurvature per- 
turbations. In general, then, at the end of the slow roll 
phase the value of the background inflaton field will have 
acquired a spatial dependence 
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where to denotes the epoch after termination of the slow- 
roll phase but before the decay process. It is necessary to 
point out that up to this point in the evolution history 
of the Universe the only difference between the single- 
field and the multi-field scenarios consists of the fact that 
field perturbations in the former case are of the adiabatic 
type only, while in the latter case they consist of both 
adiabatic and isocurvature perturbations. 

Moving on to the next phase in the evolution history of 
the Universe, it is very important to note that the quan- 
tum fluctuations imprinted onto the value of the inflaton 
fields during the slow-roll phase also represent, through 
Eq. |(2J), perturbations in the initial conditions for the 
evolution of the background fields during the decay 
phase. Similar to a chaotic mechanical system, if the 
potential V{(j)) is characterized by a global broken sym- 
metry, the system of differential equations Eqs. (|lallb() 
then exhibits a sensitive dependence on the initial condi- 
tions. This fact, in turn, guarantees that different initial 
conditions will correspond to evolution histories for the 
background fields 0(t) for t > to which do not just differ 
from one another by a simple time translation. Then to 
summarize, the broken symmetry of the potential is suffi- 
cient to guarantee that the field perturbations generated 
during the slow-roll phase will correspond to field trajec- 
tories (/)(t) that will differ substantially in the subsequent 
phase of inflaton decay. 

Moving on to consider the decay process, it is clear that 
if the inflaton fields (p have to decay into some other par- 
ticles X, then the two must necessarily be connected by 
an interaction Lagrangian C^,/,. Furthermore, this inter- 
action Lagrangian necessarily must depend on the value 
of the fields involved, (p and x- Denoting by the co- 
moving number density of x particles produced during 
the decay process, it is then trivial to note that such a 
quantity will be a (usually very complicated) functional 
of the evolution histories of the background fields. But 
since the field trajectories will depend on the initial con- 
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editions set at time to, then formally we may write 



(3) 



In a multi-field scenario the perturbations generated 
during the slow-roll phase will appear as perturbations of 
the initial conditions for the evolution of the background 
fields during the decay stage. This will in turn result in 
fluctuations in the comoving number density of particles 
produced 



(4) 



where ti > to denotes the end of the decay process. 

To determine the curvature perturbations produced 
during such a decay process, it is sufhcient to note that 
the generated energy density is proportional to the co- 
moving number density of produced particles, ~ "-x- 
Assuming for simplicity complete decay, an estimate of 
the curvature perturbation generated during this phase 
is then 



n 



Px 



(5) 



X 



where in the last step the spatially flat gauge has been 
assumed and a is a proportionality constant whose nu- 
merical value depends on the redshifting of the particle 
produced. 

The level of the non-Gaussianity produced is usu- 
ally specified using the non-linearity parameter /nl, 
which determines the non-Gaussian contribution to the 
Bardeen potential $ according to 



(6) 



where is the Gaussian part of the Bardeen potential. 
It is then possible to connect $ to C through^ 



3 ^ 3a 

5 5 Uy ' 



(7) 



where in the last step the expression for the curvature 
perturbation ^ produced in the present model has been 
used. 

It seems important to stress here that the only cru- 
cial assumption made so far concerns the global broken 
symmetry of the potential. If the symmetry is unbroken, 
it is possible to see that the perturbations of the initial 
conditions would produce evolution histories of the back- 
ground fields (j){t) which would only differ by a time trans- 
lation: the curvature perturbations thus produced would 
simply represent a gauge artifact and could be gauged 



^ It is important to note that Eq. connecting $ and is valid 
only at linear order. It is nonetheless correct to use it provided 
that the value of the non-Gaussianity parameter is larger than 
unity, which will be the case in what follows. 



away by a suitable choice of the slicing and threading. 
The only other assumption made is that the inflaton de- 
cays into the x field through the non-perturbative pro- 
cess of preheating, but this assumption is not crucial to 
the final conclusions that fluctuations in the initial con- 
ditions at the beginning of the decay stage would yield 
fluctuations in the energy density and therefore in the 
curvature. 

Curvature perturbations produced through this mech- 
anism can be characterized by significant degree of non- 
Gaussianity. The next section contains a simple example 
which shows how it is possible to determine the value of 
the non-Gaussianity parameter /nl in a specific case. 



III. AN ANALYTIC EXAMPLE: THE BROKEN 
U(l) CASE IN THE INSTANT PREHEATING 
SCENARIO 

Because of the general characteristics of the process 
being considered - a sensitive dependence of the field- 
space trajectories on the initial conditions and the non- 
perturbative nature of the preheating process - the gen- 
eral evaluation of the curvature perturbation and of the 
non-linearity parameter generated usually require a nu- 
merical approach. In this section the instant preheating 
model of Felder et al. is considered in order to obtain 
an analytic estimate of /atl- It is important to stress, 
however, that while the choice of the specific preheating 
model is in part dictated by computational convenience, 
the results of the present section can be easily general- 
ized to take into account different decay processes and 
background-field spaces. 

To proceed further it is necessary to specify two de- 
tails that so far have been left undetermined: the char- 
acteristics of the background-field space (dimensionality, 
potential, and broken symmetry) and the details of the 
preheating mechanism. The former determines the field 
trajectories (j){t), while the latter allows the specification 
of the functional relating the trajectories to the comoving 
density of particles produced, n-^[(t>{t)]. 



1. The Choice of the Background Field Landscape 

Proceeding as in Ref. 3], the dimensionality of the 
field space is chosen to be 2, so that the inflaton- field 
landscape is described by <j)i and (^2. This choice re- 
flects the fact that curvature perturbations and possibly 
non-Gaussianity can arise in the simplest multi-field in- 
flationary scenario. The complex field can then be 
introduced: 



= 01 + i02 = |0|e" 



(8) 



The potential y((^i,(^2), which is characterized by a 
slightly broken U(l) symmetry, can be modeled close to 
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its minimum by 



[l + x) 



(9) 



where x represents a measure of the symmetry break- 
ing. This choice for the potential landscape allows one 
to determine analytically the field trajectories (t){t) if the 
Hubble expansion is neglected for the short time period 
during which the decay process occurs (see the Appendix 
of Ref. 3J). 

Finally, the initial conditions for the background dy- 
namics set at the end of the slow-roll phase - and 
fluctuations thereof - can be equivalently expressed as 
[(/)i(io), 02(io)] or as [|(/)|(io), 6'(to)]- However, since the 
aim of the present work is to analyze the conversion of 
isocurvature perturbation into adiabatic perturbations 
during the decay stage, and since the form of the po- 
tential is such that the motion of the field is almost com- 
pletely in the radial direction, in what follows the latter 
choice will be adopted together with the shorthand no- 
tation [| 00 1, 6*0] • 



2. The Choice of the Inflaton Decay Process 

In this analysis, the instant preheating model of Folder 
et al. '\V^ is chosen as the process through which the in- 
flaton decays into lighter particles. Such a choice is made 
both because it is not unreasonable to think that the pre- 
heat field X is coupled to other fields into which it can 
perturbatively decay, and because it also yields a conve- 
nient analytical expression for the functional n^[0(t)]. 

The complex inflaton field </> is coupled to the preheat 
field X through the standard preheating interaction La- 
grangian C^tp — —^g\(j)\^x^- The preheat field x is also 
coupled to a massless fermion field tp by the interaction 
Lagrangian = hipil'Xj so that the perturbative de- 
cay x ~^ V'V' can take place. ^ Applying the results of 
Kofman et al. ^[ and of Felder et al. Il2l | . the comoving 
number density of x particles produced during the first 
half-oscillation of close to its minimum can be deter- 
mined to be 



3/2 



87r3 



■ exp 



(10) 



where is the instant when the inflaton reaches the 
minimum of the potential V{(f>) along its trajectory and 
1 0*1 and 1 0*1 represent the field velocity and distance 
from the absolute minimum of the potential at such an in- 
stant. From the interaction Lagrangian C^,f, it is possible 
to note that as the inflaton moves away from the bottom 



of the potential, the preheat field is endowed with an in- 
creasing effective mass = g\4>\. The x particles thus 
produced then acquire a large mass and decay into tp par- 
ticles with a decay rate T^^^^ = h?m^/S-K = h?g\(f)\/%iT 
that is proportional to |0|, which is increasing as the 
background inflaton moves away from the bottom of the 
potential. Depending on the values of the coupling con- 
stants, the whole decay chain 4> x ~^ 4' can proceed 
very efficiently, turning all the energy initially stored in 
the fields into the fermions -0 in a single half-oscillation 
of the inflaton around the minimum of its potential. The 
fermions then rapidly thermalize and complete the re- 
heating process. 

The fluctuations in the final energy density (and thus 
in the CMB temperature) can then be traced back to 
the fluctuations in the number density of preheat par- 
ticles produced during the decay process, which are in 
turn caused by the isocurvature perturbations generated 
during the slow-roll phase. 



IV. THE CALCULATION OF THE 
NON-LINEARITY PARAMETER 

To obtain an estimate of the curvature perturbation 
through Eq. jsj, it is necessary to find an expression for 
5n^/n^. From Eq. H10|l it is straightforward to get 



5n^ 



7i"5l0*P\ ^10*1 27r5|0,|" 



10* 



10* 



10*1 10* 



(11) 



Moving on to consider the background-field trajecto- 
ries, it is possible to note that the choice of the potential, 
Eq. @, allows one to derive the following approximate 
values for |0*| and |0*| as functions of the initial condi- 
tions [100 1, ^0] E^nd of the value of the symmetry breaking 
parameter x (see the Appendix of Ref. Q): 



|0.(|0o|,eo;:^)| « ^^|sin(20o)|, 



(12a) 



|04|0o|,0o;a:)| « m|0o|^l-a;sin^(^o). (12b) 

Let's then note that the form of the potential, Eq. 0, 
is such that the motion of the inflaton is almost com- 
pletely in the radial direction. Perturbations in the ra- 
dial direction therefore correspond to adiabatic pertur- 
bations, while perturbations in the azimuthal direction 
- leaving the value of the potential unaltered - corre- 
spond to isocurvature perturbations.'^ Since the aim of 
the present work is to analyze the conversion of isocurva- 
ture perturbations into adiabatic perturbations and the 



^ For a more thorough analysis of this point, involving the fact 
^ Given the large mass of the other fields involved in the decay that the field trajectory defines a local coordinate system in field 

process, the assumption of a massless i/i is actually immaterial. space, we refer the reader to the discussion given in Ref. ■ 
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resulting curvature perturbations, let's then assume that 
the slow-roll phase of inflation was such thaf* 



9o» 



(13) 



To obtain an estimate of the value of the non-Gaussianity 
parameter /nl, the Bardeen potential $ needs to be ex- 
pressed as 



Then, recalling Eq. JHl, it is possible to identify 



Inl — 



f2{0o){Seo? ^ /2(go) 

f?ido)' 



(14) 

(15a) 
(15b) 



It is also necessary to point out that while in Ref. 0] we 
performed the calculation also including terms of order 
and , in what follows only terms linear in the symmetry 
breaking parameter will be considered. Let's then notice 
that the term 



7r35x2|(/.o|sin2(20o) 



(16) 



8m\ 1 — X sin 



is of order x^, which then means that every term that 
contains it as a multiplying factor can safely be neglected 
at linear order in x. This then necessarily implies that 



$ « - 



9a I 



10 



(17) 



From Eqs. (|12all2b|l it is then possible to derive the sec- 
ond order expressions for 



10* 



cos(2^ 



9o^ 



(18a) 



x 



10*1 2[1 
cos(26'o) + 



— X sin^(0o)] 

j:sin2(26'o) 
4[l-j;sin2(6'o)] 



sin(26'o) 



(18b) 



where the expression for (5|0, |/|0, | has been added for 
completeness. We then have 



$ = 



9ax 



20 [1 - 
cos(26io) 



xsin^(6'o)] 

xsin^(2^^o 



sin(26lo)(56l 



4[1 



X sm' 



(19) 
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FIG. 1: Value of /nl for values of the symmetry breaking 
parameter of x = 10"^ (solid curve) and x — 10~^ (dashed 
curve). It seems important to stress that only a very small 
interval in the initial conditions 6o will yield insignificant non- 
Gaussianity. 



which then leads to the identification 
9ax 



f2iOo) = 



20 [l- xsin^ieo)] 
9ax 



sin(26'o), 



(20a) 



20 [l- j:sin2(6lo)] 

foa \ , 2:sin^(2e'o 
cos(20o) 



4[l-a;sin2(6'o)] 
It is then immediate to find that 
20 [1 -a;sin2(6'o)] 



(20b) 



Inl — 



9a a;sin^(26'o) 

^ X sin2(26'o) 
cos(26lo) + -- ^ ^ 



4 [1- j;sin2(6'o)] 



20 cos(26lo) 
9aa;sin2(26lo)' 



(21) 



where in the last expression only the leading term has 
been considered.^ It is then very interesting to note this 
particular model is characterized by a value of the non- 
linearity parameter which is inversely proportional to the 
value of the symmetry breaking: in general slightly bro- 
ken symmetries should then produce mild curvature per- 
turbations characterized by large non-Gaussianity. 



* This happens if the potential in the slow-roll region is such that 
the effective mass for excitation of quantum fluctuations of the 
inflaton in the direction of motion is larger than the Hubble 
parameter during slow-roll. 



It is in fact important to note that if we wanted to perform a 
calculation exact up to a;" order we should also consider a term 
arising from the (5|0* 1/10, | factor. 
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V. GENERAL RESULT FOR /nl AND 
CONCLUSIONS 



characterized by a global broken symmetry. In particu- 
lar, it is straightforward to obtain the following condition 



The results of the previous section were obtained un- 
der very specific assumptions: a broken C/(l) symmetry 
and the instant preheating mechanism. It is nonetheless 
possible to generalize the previous approach to obtain a 
general formula to estimate the value of the non-linearity 
parameter of the curvature perturbations produced in a 
multi-field scenario by a generic decay process. 

As it was shown in Ref. 0], the initial conditions per- 
turbation S(f)(tQ,x) can in fact be decomposed into an 
adiabatic component i5(/)|| parallel to the field velocity 

M|| = d(f>/dt and into an isocurvature component S(f>± per- 
pendicular to it and directed along the unit vector uj_. 
As in Ref. Isl, it is assumed here as well that the adia- 
batic perturbations generated during the slow-roll phase 
of inflation are actually negligible compared to the isocur- 
vature component 



> 



(22) 



This immediately leads to the following estimate for the 
Gaussian curvature perturbation 



3a dn^ 
5nv 06 I 



(23) 



It is important to stress that the functional relating the 
comoving number density of particles produced during 
the decay process to the background trajectory is in gen- 
eral very complicated, usually requiring numerical solu- 
tions to formulate quantitative predictions. Furthermore, 
if the potential is characterized by a broken symmetry 
then the background field trajectories may also be char- 
acterized by a sensitive dependence on the initial con- 
ditions and therefore may exhibit chaotic behavior. All 
these factors then conspire to make the evaluation of $ 
rather involved. 

Recalhng Eqs. (OEJ) and expanding 5n^/n^ to second 
order and neglecting the terms containing powers of 5(j)\\ , 
the Bardeen potential turns out to be 



3a 



(24) 



which then allows one to obtain the following general 
expression for the non-linearity parameter 



NL 



d'^n^/d4>']^ 



(25) 



It is then clear that depending on the relative mag- 
nitudes of the derivatives of n^^, significant or negligible 
non-linearity can be produced by a multi-field potential 



I d'^Uy ^ 2a f dn 



d(f)\ 



(26) 



which, given the chaotic nature of the background field 
trajectories and the exponential production of particles 
in the resonances characteristics of the preheating pro- 
cesses, can possibly be satisfied by a large class of mod- 
els. 

A few comments are finally in order. The level of non- 
Gaussianity predicted by the model at hand can be rather 
sizable, values oi Jnl ~ 100 can be easily achieved. 
As we have already mentioned in the Introduction, in 
the computation of the nonlinearity parameter /^vl we 
have neglected the nonlinear gravitational contribution 
coming from the post-inflationary evolution. Since they 
amount to a contribution to Jnl of order unity 0, it 
follows that the estimate of /jv l will be valid only in the 
regime f^L ^1- In other words, the validity of the 
estimate requires that the global symmetry is slightly 
broken, that is x ^ 1. The fact that the level of non- 
Gaussianity is inversely proportional to the parameter x, 
which gives a measure of the breaking of the global sym- 
metry within the infiaton sector, is reminiscent of what 
happens for the curvaton scenario [l^ . There the non- 
Gaussianity is inversely proportional to the parameter r 
which measures the energy density contribution of the 
curvaton to the total one at the time of curvaton decay. 
Since present WMAP data impose a bound I/nl] ^ 100, 
this translates into a weak bound on the breaking param- 
eter 



> 10-2 tan (26lo) sin {260) . 



(27) 



Since future tests of the angular bispectrum of tem- 
perature anisotropy alone could detect non-Gaussianity 
as small as I/nl] ~ 5 0, it would be interesting 
to generalize the computation performed in this pa- 
per to include the contribution of gravitational nonlin- 
earities which render the nonlinearity parameter 
momentum-dependent. As explicitly shown in Ref. 
this dependence augments the non-Gaussian signal at 
large multipoles and may be therefore crucial in testing 
large symmetry breaking patterns in the inflaton sector. 
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